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ABSTRACT 

A  Study  of  the  performance  parameters  for  a  journal  bearing  of  finite 
length  lubricated  with  micropolar  fluids  is  undertaken.  Results  indicate  that 
a  significantly  higher  load  carrying  capacity  than  Newtonian  fluids  may  result 
depending  on  the  size  of  material  characteristic  length  and  the  coupling 
number.  It  is  also  shown  that  although  the  frictional  force  associated  with 
micropolar  fluid  is  in  general  higher  than  that  of  a  Newtonian  fluid,  the 
friction  coefficient  of  micropolar  fluids  tends  to  be  lower  than  that  of  the 
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nondimensional  film  thickness  =  - 
bearing  length,  m 


nondlmensional  characteristic  length  of  micropolar  fluid  =  ^ 
coupling  number 
hydrodynamic  pressure,  N/m^ 

nondimensional  pressure  =  ^  jj 

leakage  flow  rate,  m3/s 

nondimensional  leakage  =  r  Q 

cUR 

journal  radius,  m 

component  of  load  along  the  r  ■'ial,  N 


nondimensional  load  component  =  - ; 

pULR' 


component  of  load  along  the  tangential  direction,  N 

W  c2 

nondimensional  load  component  =  — 1 — « 

pULR4 

resultant  of  load,  N 

Wc2 

nondimensional  load  =  - ~ 

yiULR4 


coordinate  system,  m 

y 

nondimensional  coordinates  =  £ 
micropolar  viscosity  coefficient,  N-s 
eccentricity  ratio,  c  =  e/c 
nondimensional  coordinate,  x/R 
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0Cav  angular  coordinate  where  the  film  cavitates 
k  micropolar  viscosity  coefficient,  N  s/m2 

A  characteristic  length  of  micropolar  fluids,  m 

p  viscosity  of  the  Newtonian  fluid,  N  s/m2 

<t>  function  defined  in  Eq.  (4) 

$  normalized  form  of  4> 

<p  attitude  angle,  rad 

INTRODUCTION 

In  recent  years  there  has  been  a  considerable  amount  of  interest  in  the 
lubricating  effectiveness  of  non-Newtonian  fluids.  Generally  speaking,  in 
non-Newtonian  lubricants  the  shear  stresses  do  not  always  remain  proportional 
to  the  shear  rates,  particularly  for  high  shear  rates.  To  model  such  fluids, 
various  mathematical  formulations  such  as  the  power  law  and  Ree-Ering  models 
have  been  proposed.  This  paper  presents  a  study  of  the  non-Newtonian  behavior 
of  the  lubricating  oil  (in  a  finite  length  journal  bearing)  from  the  modern 
continuum  point  of  view.  In  continuum  mechanics,  the  continuous  media  is 
regarded  as  a  set  of  structured  particles  possessing  individual  mass  and 
velocity.  These  continuum  particles  are  assigned  a  sub-structure  which  can 
translate,  rotate,  or  even  deform  independently.  In  other  words,  the 
microvolume  is  considered  to  have  independent  microvolume  elements  that  possess 
individual i ty. 

Theoretical  work  dealing  with  the  mechanics  of  fluids  with  microstructure 
began  with  the  work  of  Jeffery  (Ref.  1)  in  1922  who  studied  the  problem  of  a 
Newtonian  fluid  motion  in  the  presence  of  suspended  ellipsoids.  His  results 
indicated  that  the  viscosity  of  the  fluid  with  such  particles  was  higher  than 


,  **»  4  I  i- 1  |*t  l1.  O.  in  aV if.tiUi<>,lttM.*>l.llt4.,»«  nlMALM  M.*! 


that  of  its  base  Newtonian  fluid.  Later  in  1957,  Prager  (Ref.  2)  derived 
constitutive  equations  for  stress  in  a  noninteracting  (dilute)  dumbbell 
particle  suspension.  Since  unlike  the  classical  fluid  mechanics  the 
microcontinua  is  regarded  as  sets  of  structured  particles,  the  theory  required 
a  complete  re-examination  of  the  classical  concepts  that  has  to  be  augmented 
with  additional  balance  laws  and  constitutive  relations. 

The  general  theory  of  microcontinua  is  attributed  to  Eringen  who  developed 
the  fluid  mechanics  of  deformable  microelements  (Ref.  3)  in  an  article  entitled 
"Simple  Microfluids."  The  Eringen  definition  of  a  simple  microfluid  is  a 
viscous  medium  whose  behavior  and  properties  are  affected  by  the  local  motion 
of  particles  in  its  microvolume.  These  fluids  are  characterized  by  22 
viscosity  and  material  constants  and  when  applied  to  flow  problems  the  result 
is  a  system  of  19  partial  differential  equations  with  19  unknowns  that  may  not 
be  amenable  to  solutions.  Eringen  (Ref.  4)  subsequently  introduced  a  subclass 
of  fluids  which  he  named  micropolar  fluids  that  ignores  the  deformation  of  the 
microelements  but  still  allows  for  the  particle  micromotion  to  take  place. 

Since  then  many  other  papers  have  appeared  in  the  literature  that  present 
different  views  of  the  subject  based  on  the  continuum  mechanics.  A 
comprehensive  review  of  these  theories  is  provided  by  Ariman  et  al.  (Ref.  5). 

As  for  the  application  of  the  microcontinuum  theory  to  real  flows,  the 
blood  rheology  appears  to  be  a  very  promising  area  (Refs.  6  and  7).  Continuum 
approach  to  blood  flow  in  the  arterioles  takes  the  rotation  and  deformation  of 
red  blood  cells  into  account.  Such  consideration  may  be  particularly  important 
because  of  the  relative  size  of  the  cells  as  compared  with  the  diameter  of  the 
arterioles.  Another  area  where  microcontinuum  theory  may  be  applied  is  in  the 
study  of  rheological  behavior  of  complex  microstructures  such  as  liquid 
crystals  which  are  known  to  possess  complex  molecular  structure  and  orientation 
patterns  (Ref.  8). 
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Microcontinuum  theory  In  general  may  also  be  applied  to  the  study  of  the 

non-Newtonian  behavior  of  certain  polymer  solutions  such  as  polyisobutylene 

which  are  known  to  have  long  chains  of  molecular  structure.  To  study  such 

fluids.  Hand  (Ref.  9)  put  forward  a  theory  of  so  called  anisotropic  fluids  and 

presented  theoretical  results  with  comparison  to  those  measured  experimentally 
* 

Hand  showed  that  theories  due  to  Jeffery  (Ref.  1)  and  Prager  (Ref.  2)  are 
contained  in  his  anisotropic  fluids  as  special  cases  and  that  the  theory  can 
be  applied  to  fluids  with  microstructure  such  as  polyisobutylene  that  are 
spherical  at  rest  but  elongate  under  shear  field  as  do  immiscible  drops 
suspended  in  fluids.  Comparing  results  with  experimental  measurements.  Hand 
reports  good  agreement  for  the  dependence  of  apparent  viscosity  and  normal 
stress  particularly  at  low  shear  rates.  His  success  is  indicative  of  the 
potential  of  applicability  of  modern  continuum  theory  for  non-Newtonian  fluids 

Theoretical  works  on  the  non-Newtonian  bearings  from  microcontinuum 
viewpoint  began  with  the  work  of  Allen  and  Kline  (Ref.  10).  They  stated  that 
in  addition  to  polymeric  fluids  with  long  molecular  chains,  the  continuum 
model  can  be  used,  as  a  first  approximation,  for  modeling  bearings  that  are 
contaminated  with  dirt  and  particles.  Under  such  conditions,  the  lubricant 
can  be  considered  as  a  fluid  suspension.  They  studied  the  one-dimensional 
slider  bearing  problem  and  concluded  that  the  theory  predicts  a  higher  load 
capacity  and  a  lower  friction  coefficient  compared  to  those  of  a  Newtonian 
fluid,  thereby  signifying  an  improvement  in  the  lubrication  of  a 
one-dimensional  bearing  lubricant  with  micropolar  fluids. 

Since  that  article  was  published,  many  other  papers  appeared  in  the 
literature  that  dealt  with  the  lubrication  of  one-dimensional  bearings.  To 
put  things  in  perspective,  the  following  references  can  be  mentioned:  for 
slider  or  step  bearings  (Refs.  11  to  17),  thrust  bearings  (Refs.  18  and  19), 


journal  bearings  (Refs.  20  to  25),  squeeze  film  bearings  (Refs.  26  to  30), 
rolling  element  and  elastohydrodynami c  lubrication  (Refs.  31  and  32). 


Governing  Equations 

The  general  form  of  the  governing  equations  for  the  steady  state  motion 
of  micropolar  fluids  can  be  written  in  tensor  notation  as  follows  (Ref.  31). 
conservation  of  mass: 

7  •  (pV)  =0  (1) 

Conservation  of  linear  momentum: 

(X  +  2|i)  77  •  7  -  ^  (2p  +  tc)7x7  +  KVxu-7ir  +  pFg  =  0  (2) 

Conservation  of  angular  momentum: 

(a  +  3  +  y)V7  •u-y7x7xu+»c7x\/-  2kv  +  p  £q  =  0  (3) 

where  V  i s  the  velocity  vector  and  u  is  the  mi crorotational 

velocity  vector,  p  is  the  mass  density,  ir  is  the  thermodynamic  pressure,  b 

Fb  Is  the  body  force  per  unit  mass,  £q  is  the  body  couple  (moment)  per  unit 
mass,  \i  and  X  are  the  familiar  viscosity  coefficients  of  the  classical 
fluid  mechanics,  while  a,|3,y,  and  k  a-e  additional  viscosity  coefficients 
for  micropolar  fluids. 

Equation  (2)  has  the  familiar  form  of  the  Navier-Stokes  equation  but  it 

is  coupled  with  Eq.  (3)  which  essentially  describes  the  motion  of  the  particles 

inside  the  microvolume  as  they  undergo  mi crorotational  effects  represented  by 
the  microrotational  vector  u.  For  fluids  with  no  microstructure,  this 
parameter  vanishes.  For  Newtonian  fluids,  Eqs.  (2)  and  (3)  decouple  since 
•c  *  0.  These  equations  can  be  simplified  considerably  for  incompressible 
fluids  (p  =  constant)  in  the  absence  of  body  forces  and  inertial  forces  since 


7  •  V  =  0  and  Fg  and  £g  are  assumed  to  be  negligible.  Furthermore,  the 
thermodynamic  pressure  ir  can  be  replaced  by  P  since: 


~  9E  ' 

.ap_1J 


where  E  is  the  internal  energy. 


Having  made*  the  aformentioned  assumptions,  a  so-called  generalized 
Reynolds  equation  for  micropolar  fluids  can  be  derived,  see  Ref.  14.  The 
resulting  equation  is  given  below  (Ref.  25): 
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3x  .y 


fT  A/M  a  *^3P  8_  hi  A  h,  3£  U  3h 
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The  parameter  y  is  the  viscosity  of  the  base  fluid  as  in  the  case  of 
Newtonian  fluids  and  k  and  y  are  two  additional  viscosity  coefficients  for 
micropolar  fluids.  These  viscosity  parameters  are  grouped  in  the  form  of  two 
parameters  N  and  A.  N  is  a  dimensionless  parameter  called  the  coupling 
number,  for  it  characterizes  the  coupling  of  the  linear  and  angular  momentum 
equations.  When  N  is  identically  zero,  the  equations  of  linear  and  angular 
momentum  are  decoupled  and  the  equation  of  the  linear  momentum  reduces  to  the 
classical  Navier-Stok.es  equation.  The  parameter  A  is  called  the 
characteristic  length  for  it  characterizes  the  interaction  between  the 
micropolar  fluid  and  the  film  gap.  This  parameter,  carrying  the  dimension  of 


length,  is  a  function  of  the  size  of  the  lubricant  molecule.  As  A  approaches 
zero,  the  effect  of  microstructure  becomes  less  important  and  in  the  limiting 
case,  when  A  vanishes,  the  function  <t>  becomes  identical  to  1/12  and  the 
Reynolds  equation  reduces  to  its  classical  form. 


Nondimensional ization 


a=£  e  =  *  n  Y  R_h  p_£Rc 

m  A  y  Ry‘L  n‘c  K  “  p  U  R 

Then  the  generalized  Reynolds  equation  for  micropolar  fluids  becomes: 


3_  If3  3p|  /R\2  3  3  5  3P]  _  ,  3fi_ 

ae  L  *  aej  +Vi/  ay  Lfi  ay  J  =  °  ae 


where : 


*  -  I  +  12  -  6N 

jtr2 

m 


Boundary  conditions  are: 

P  =  0  at  e  =  0 

AS 

P  -  55  ■  0  at  e  -  ecav  (6> 

P  =  0  at  y  =  0  and  y  =  1 

where  9  =  9cav  is  the  angular  coordinate  at  which  the  film  cavitates. 

Equation  (5)  is  a  fully  elliptic  partial  differential  equation  for  which 
no  analytical  solution  exists.  Analytical  solutions  for  a  number  of  bearing 
configurations  are  available  as  cited  earlier  for  either  infinitely  long 
bearings  (L/D  >>  1)  or  short  bearings  (I/O  <<  1).  To  evaluate  the  performance 
of  a  bearing  for  a  finite  length,  one  must  resort  to  numerical  schemes  for 
solution.  Equation  (5)  was  solved  using  the  finite  difference  method:  the 
equation  was  discretized  and  the  difference  equations  were  solved  on  a  computer 
using  the  Successive-Over-Relaxation  method  (S.O.R.).  Convergence  was  assumed 


when  the  relative  maximum  error  between  two  successive  iterations  fell  below  a 
tolerance  value  of  10-3.  The  Reynolds  cavitation  boundary  conditions  were 
satisfied  by  setting  all  negative  pressures  equal  to  zero  during  the 
computations.  Once  the  pressure  distribution  was  evaluated,  the  following 
performance  parameters  were  computed  numerically. 

Load  capacity  and  attitude  angle: 


WRc2  f  p6cav 

Wp  =  - =•  =  J  J  p  cos0  d0  dy 

pLILR^  0  0 

W  c2  1  0cav 

R  =  -*—5  =  1 1  p  sin©  d0  dy 

*  pULRz  0  0 


(7) 


(8) 


h.-SsL. 

pULR^ 


Leakage  flow  rate: 


Friction  force: 


where: 


<p  =  tan 


o  =  Q  =  Jecav  *  fi3  |E  de 

CUR2  0  9y 


F  =  -~F  =  J  A  d0  +  J  A  {—■  d6 

^  0  9cav  ^ 


(9) 


(10) 


(11) 


(12) 


2  d0  +  r  2N f cosh ( NLH )  -  1] 
L  L  sinh(NLPi)  J 

Friction  coefficient: 


(13) 


Equations  (7)  to  (13)  were  numerically  evaluated  using  a  Simpson's  integration 


scheme . 


RESULTS 


To  study  the  lubricating  effectiveness  of  fluids  that  exhibit 
micropolarity,  a  parametric  study  was  performed.  Results  are  obtained  as  a 

function  of  two-nondimensional  parameters,  N  and  fi.m  which  represent  the 

* 

characteristics  of  micropolar  fluid.  Unless  otherwise  stated,  the  computations 
pertain  to  a  finite  journal  bearing  of  L/D  =  1. 

Pressure  Distribution  and  Load  Capacity 

Figure  1  presents  a  number  of  three-dimensional  pressure  profiles  for  a 
steadily  loaded  journal  bearing  at  an  eccentricity  ratio  of  e  =  0.5. 

Figure  1(a)  depicts  the  Newtonian  case  while  Fig.  Kb)  and  (c)  show  the  effect 
of  micropolar  lubricant  on  the  pressure  distribution.  In  Fig.  Kb),  the 
nondimensional  characteristic  length,  dm,  was  held  constant  at  dm  =  9.0  while 
the  coupling  number  parameter,  N2 ,  was  varied  from  0.1  to  0.9.  Computations 
indicate  that  although  the  shape  of  the  pressure  distribution  remains  the  same, 
the  magnitude  of  the  pressure  and  thus  the  load  capacity  varies  significantly 
as  a  function  of  the  coupling  number.  From  thermodynamic  considerations 
parameter  N  is  limited  to  0  <  N  <  1,  see  Eringen  (Ref.  3).  At  high  values 
of  N  the  effect  of  microstructures  becomes  significant  whereas  at  low  values 
of  N  the  individuality  of  the  substructure  is  much  less  pronounced.  In  the 
limiting  case,  as  N  approaches  zero,  the  micropolarity  is  lost  and  the 
lubricant  is  considered  to  behave  as  Newtonian.  Under  such  a  limiting  case, 
the  equation  of  linear  momentum  reduces  to  that  of  the  Navier-Stokes . 

The  effect  of  the  nondimensional  characteristic  length  on  the  pressure 
distribution  is  depicted  in  Fig.  1(c).  As  shown,  at  low  dm,  the  effect  of 
microstructure  is  very  pronounced.  Physically,  lower  values  of  dm  correspond 


to  larger  values  of  the  characteristic  length  (Q.m  =  c/A,  where  c  -  the  bearing 
clearance  -  is  assumed  to  remain  constant).  Small  &m,  therefore,  indicates 
that  characteri sti c  length  of  the  substructure  is  large  compared  to  that  of  the 
clearance  dimension.  At  high  values  of  &m,  the  effect  of  substructure  becomes 
less  and  less  significant  and  in  the  limiting  case  when  Q.m  approaches 
infinity,  the  effect  of  individuality  of  the  microstructure  is  lost  so  that  the 
pressure  distribution  approaches  that  of  the  Newtonian. 

Variation  of  the  load  carrying  capacity  as  a  function  of  the 
eccentricity  ratio,  e,  is  depicted  in  Fig.  2  for  various  values  of  N,  keeping 
the  nondimensional  characteristic  length  constant,  2.m  =  9,  for  all  cases.  The 
enhancement  in  the  load  capacity  for  higher  values  of  the  coupling  number  is 
pronounced.  A  possible  explanation  for  such  phenomena  can  be  attributed  to  the 
increase  in  the  viscosity  due  to  the  additives,  characterized  by  micropolar 
fluids.  There  are  some  experimental  works  that  support  this  view.  Henniker's 
extensive  review  (Ref.  33)  entitled,  "The  Depth  of  the  Surface  Zone  of  a 
Liquid,"  (Ref.  33)  attest  to  the  evidence  of  abnormally  high  viscosity  in  the 
vicinity  of  a  solid  surface.  He  refers  to  some  experimental  work  by  Deryagin 
et  al.  (Ref.  34)  who  found  the  viscosity  of  a  turbine  oil  within  5000  A  of 
solid  wall  to  increase  tenfold  when  aluminum  naphthionate  of  up  to  2  percent 
were  added  in. 

Figure  3  depicts  how  the  load  carrying  capacity  varies  as  a  function 
of  2.m  for  a  number  of  ,  keeping  the  eccentricity  constant.  It  clearly 
shows  an  appreciable  rise  in  the  load  capacity  at  small  values  off  Q.m  and 
that  the  load  approaches  to  that  of  the  Newtonian  fluids  in  the  limiting  case 
as  -»  ®.  These  results  are  consistent  with  the  analytical  solutions  due  to 
Cowin  (Ref.  35)  who  showed  that  for  polar  fluids  in  couette  flow  and  Poiseuille 
flow  the  effective  viscosity  ue  approaches  the  Newtonian  fluid  value  p  as 
*  00  ■ 


Friction  Coefficient 


Variation  of  the  coefficient  of  friction  as  a  function  of  eccentricity  is 
shown  in  Fig.  4  for  a  number  of  N2  while  Q.m  =  9.0.  Also  shown  is  the 
corresponding  variation  in  friction  coefficient  for  the  Newtonian  fluid  case. 
This  figure  indicates  that,  under  the  conditions  simulated,  micropolar  fluids 
exhibit  a  beneficial  effect  in  that  the  coefficient  of  friction  is  lower  than 
the  Newtonian  fluids.  Furthermore,  as  the  coupling  number  increases,  a  lower 
friction  coefficient  may  result.  To  further  explore  this  effect,  the 
coefficient  of  friction  was  plotted  as  a  function  of  the  for  a  number  of 
N2  as  presented  in  Fig.  5.  This  figure  clearly  shows  that  micropolar  friction 
coefficient  approaches  that  of  Newtonian  results  as  Hm  ■*  0  and  as  Q.m  ■*  °°  or 
as  N2  ->  0.  The  most  interesting  result  is  that  there  exists  an  optimum 
material  length  at  which  the  friction  coefficient  is  minimum. 

It  should  be  noted  that  although  the  friction  coefficient  of  lubricant 
micropolarity  tends  to  be  lower  than  that  of  the  Newtonian  fluids,  the  friction 
force  tends  to  be  higher  as  depicted  in  Fig.  6.  The  reason  for  this  is  of 
course  due  to  the  load  capacity  of  the  micropolar  fluids  increasing  more  than 
the  corresponding  increase  in  (Fig.  3)  which  more  than  compensates  for  the 
frictional  force.  Similar  behavior  were  predicted  by  analytical  work  of  Allen 
and  Kline  (Ref.  10)  for  a  one-dimensional  slider  bearing. 

CONCLUDING  REMARKS 


A  study  of  the  lubricating  effectiveness  of  micropolar  fluids  in  a 
finite  journal  bearing  is  presented.  Micropolar  fluids  are  a  subclass  of 
microfluids  in  which  the  fluids  exhibit  microrotational  effects  i.e.,  the  fluid 
is  considered  to  possess  microstructure  that  can  rotate  independently  of  the 


microvolume  but  the  deformation  (microstretch)  is  not  allowed.  These  fluids 
have  potential  in  describing  the  effect  of  polymeric  additives  as  the  classical 
Navier-Stokes  theory  has  no  provision  for  the  effects  of  microstructure  in 
fluids. 

The  results  presented  in  this  paper  indicate  that  for  a  steadily  loaded 
finite  journal  bearing,  the  micropolar  fluids  do  indeed  exhibit  a  beneficial 
effect  in  that  the  load  carrying  capacity  is  significantly  increased  and  the 
friction  coefficient  is  less  than  that  of  the  Newtonian  lubricant. 
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